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The distribution function of the free energy fluctuations in one-dimensional directed polymers 
with (^-correlated random potential is studied by mapping the replicated problem to a many body 
quantum boson system with attractive interactions. Performing the summation over the entire 
spectrum of excited states the problem is reduced to the Fredholm determinant with the Airy kernel 
which is known to yield the Tracy- Widom distribution. 
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I. INTRODUCTION 



Directed polymers in a quenched random potential have been the subject of intense investigations during the past 
two decades (see e.g. [l|). In the most simple one-dimensional case we deal with an elastic string directed along 
the r-axis within an interval [Q,£]. Randomness enters the problem through a disorder potential V[4>(t),t], which 
competes against the elastic energy. The problem is defined by the Hamiltonian 



H 



r),V] = J dr{i[5 T 0(r)] 2 + y[0(r),r]}; (1) 

where in the simplest case the disorder potential V[(f>, t] is Gaussian distributed with a zero mean V{4>, r) = and 
the ^-correlations: 



V{4>, r)V((j)', t') = u6(t - T')8{(j) - <j)') (2) 



— i . Here the parameter u describes the strength of the disorder. Historically, the problem of central interest was the 
!> scaling behavior of the polymer mean squared displacement which in the thermodynamic limit (L — > oo) is believed 
to have a universal scaling form (</> 2 }(L) oc L?^ (where (...) and (...) denote the thermal and the disorder averages), 
with £ = 2/3, the so-called wandering exponent. More general problem for all directed polymer systems of the 
type, Eq.([T]), is the statistical properties of their free energy fluctuations. Besides the usual extensive (linear in L) 
self- averaging part foL (where /o is the linear free energy density), the total free energy F of such systems contains 
disorder dependent fluctuating contribution F, which is characterized by non-trivial scaling in L. It is generally 
believed that in the limit of large L the typical value of the free energy fluctuations scales with L as F oc L 1 / 3 (see 
e.g. In other words, in the limit of large L the total (random) free energy of the system can be represented as 

F = f L + cL 1 / 3 f (3) 

, where c is the parameter, which depends on the temperature and the strength of disorder, and / is the random 
quantity which in the thermodynamic limit L — > oo is described by a non-trivial universal distribution function P*(/). 
The derivation of this function for the system with (^-correlated random potential, Eqs.([T])-([2]) is the central issue of 
the present work. 

For the string with the zero boundary conditions at r = and at r = L the partition function of a given sample is 

Z[V] = / T>[<f>(r)] e-P H[ *' v] (4) 

J0(O)=O 

where (3 denotes the inverse temperature. On the other hand, the partition function is related to the total free energy 
F[V] via Z[V] = exp(— /3F[V]). The free energy F[V] is defined for a specific realization of the random potential 
V and thus represent a random variable. Taking the iV-th power of both sides of this relation and performing the 
averaging over the random potential V we obtain 



Z N [V] = Z[N, L] = exp(-f3NF[V}) (5) 
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where the quantity in the lhs of the above equation is called the replica partition function. Substituting here F 
f L + cL 1 ? 3 /, and redefining Z[N, L] = Z[N, L] exp{-f3Nf L} we get 



Z[N,L] = exp(-AAT/) (6) 

where A = ficL 1 / 3 . The averaging in the rhs of the above equation can be represented in terms of the distribution 
function Pl{J) (which depends on the system size L). In this way we arrive to the following general relation between 
the replica partition function Z[N, L] and the distribution function of the free energy fluctuations 

r+oo 

z[n,l] = / df p L (f) e- XN f (7) 



Of course, the most interesting object is the thermodynamic limit distribution function P*(f) — liniL-^oo Pl{/) which 
is expected to be the universal quantity. The above equation is the bilateral Laplace transform of the function Pl(/), 
and at least formally it allows to restore this function via inverse Laplace transform of the replica partition function 
Z[N,L]. In order to do so one has to compute Z[N,L] for an arbitrary integer N and then perform analytical 
continuation of this function from integer to arbitrary complex values of N. In Kardar's original solution [f|, after 
mapping the replicated problem to interacting quantum bosons, one arrives at the replica partition function for 
positive integer parameters N > 1. Assuming a large L — > oo limit, one is tempted to approximate the result by 
the ground state contribution only, as for any N > 1 the contributions of excited states are exponentially small for 
L — > oo. However, in the analytic continuation for arbitrary complex N the contributions which are exponentially 
small at positive integer N > 1 can become essential in the region N — > 0, which defines the function P(f) (in other 
word, the problem is that the two limits L — > oo and N — > do not commute [1, 0|)- Thus, it is the neglection of the 
excited states which is the origin of non-physical nature of the obtained solution. 

In the my recent paper [8j an attempt has been made to derive the free energy distribution function via the 
calculation of the replica partition function Z[N, L] in terms of the Bethe-Ansatz solution for quantum bosons with 
attractive (5-interactions which involved the summation over the entire spectrum of exited states. Unfortunately, the 
attempt has failed because on one hand, the calculations contained a kind of a hidden " uncontrolled approximation" , 
and on the other hand, the analytic continuation of obtained Z(N, L) was found to be ambiguous. 

It turns out that it is possible to bypass the problem of the analytic continuation if instead of the distribution 
function itself one would study its integral representation 

/•OO 

W(x) = / dfP*(f) (8) 



which gives the probability to find the fluctuation / bigger that a given value x. Formally the function W{x) can be 
defined as follows: 

°° f— 1) N 

W(x) = Van N \ exp(A7Vx) Z N = Vm exp[- exp(A(x - /))] = 6{f-x) (9) 

On the other hand, in terms of the replica approach the function W(x) is given by the series 

00 (—l) N 

W(x) = lim V ^ ' exp(AiVx) Z\N, L] (10) 

A^oo * — ' jv! 
N=0 

In the present paper the replica partition function Z[N, L] will be calculated (again) by mapping the replicated 
problem to the TV-particle quantum boson system with attractive interactions. Performing the summation over the 
entire spectrum of excited states the summation of the series, eq. (fTtlj) . is reduced to the Fredholm determinant with 
the so called Airy kernel which is known to yield the Tracy- Widom distribution. Originally this distribution has been 
derived in the context of the statistical properties of the Gaussian Unitary Ensemble Q while at present it is well 
established to describe the statistics of fluctuations in various random systems [Iol - fl5| which are widely believed to 
belong to the same universality class as the present model [l6l - [l8| . While this manuscript was in course of preparation 
I have learned that the exactly the same result for the system considered in this paper has been independently derived 
by P.Calabrese, P.Le Doussal and A. Rosso (l9j . 

Performing simple Gaussian average over the random potential, eq.([2|), for the replica partition function, Eq.([5|), 
we obtain the standard expression 

Z(N,L) = H VMr)e~ pHlfW (11) 

a=l J 'M0)=0 
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where 

1 C L ( N N \ 

H N [<f>] = - dr V[9A(r)] 2 -/3«V*[^(r)-^(r)] (12) 

ZJ ° \a=l a^b J 

is the TV-component scalar field replica Hamiltonian and tfi = {<fri, . . . ,<Pn}. 

According to the above definition this partition function describe the statistics of N (5-interacting (attracting) 
trajectories 4> a { T ) a U starting (at r = 0) and ending (at t = L) at zero. In order to map the problem to one- 
dimensional quantum bosons, let us introduce more general object 

*(x;i) = TT/ VMr)e~ mN W (13) 

Q=1 J0 a (O)=O 

which describes N trajectories <f> a ( T ) an starting at zero (<fi a (0) = 0), but ending at r = t in arbitrary given points 
{x\, a; at}. One can easily show that instead of using the path integral, ^(x; t) may be obtained as the solution of 
the linear differential equation 

N N 

9 t *(x;t) = —^^(x;*) + ^E 6 ( x * ~ X "M^ *) ( 14 ) 

" o=l a=£b 

with the initial condition ^(x; 0) = H^ =1 5(x a )- This is nothing else but the imaginary-time Schrodinger equation which 
describes N bose-particles of mass /3 interacting via the attractive two-body potential — 2 u8{x). The original replica 
partition function, Eq. (fTTj) . then is obtained via a particular choice of the final-point coordinates, Z(N,L) = ^(0:L) 
The spectrum and some properties of the eigenfunctions for attractive one-dimensional quantum bosons have been 
derived by McGuire [2(| and by Yang [2l[ (see also Ref. [I^,[23|]). A generic eigenstate of this system consists of M (1 < 
M < N) "clusters" {f2 Q } (a = 1, M) of bound particles. Each cluster is characterized by the momentum q a of its 
center of mass motion, and by the number n a of particles contained in it (such that Y] . n a = N). Correspondingly, 

the eigenfunction Vpqjn (xi, ...,xpf) of such state is characterized by M continuous parameters q = (qx, qu) and M 

integer parameters n = (m «m) (detailed structure and the properties of these wave functions are described in 

Ref.Q). The energy spectrum of this state is 

= ^E E^) 2 = ^E - ^E(»«-"-) ( 15 ) 

^ a = l r=l 1 a=l ^ a=l 

As the wave functions ^^/(x) can be proved to constitute the complete and othonormal set, the replica partition 
function of the original directed polymer problem can be represented in the form of their linear combination: 



N / M +oc , oo \ / M 

*(x,t) = E n / E ) ^Wt r (») e-*"i«*» 5(E^, 

M=l \q = 1 J - 00 Z7T n a = J \ a =l 



(16) 



where m) is the Kronecker symbol (which allows to extend the summation over n Q 's to infinity). Using the 

explicit form of the wave functions "J^? (x) the above expression (after somewhat painful algebra) reduces to (see Q 
for details) 

Z(N,L) =e- pNLfo Z{N,\) (17) 
where /o = ^g(3 4 u 2 — w ln(/3 3 u) is the linear (selfaveraging) free energy density, and 

Z(N,X) = Nl[ M{y (18) 



4ttA7V 



,„V 1 n V / f °° dy a dp a 2\a*n a y a TT | A V"a _ 'W _ 'W* c/v^ /v 

+ ™ ! e in n e / / 4^**.+«je j n p^^-rf w 



M=2 L a=ln c «=l*' -x 



A(n Q -np)- i(p a -pp)\ 



where Ai(t) is the Airy function, and instead of the system length L we have introduced a new parameter A = 



u L) ' . The first term in the above expression is the contribution of the ground state (M = 1), and the next 



terms (M > 2) are the contributions of the rest of the energy spectrum. 
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Using the Cauchy double alternant identity 

[ILK ~ap)(b a -bp 



= (-l) N ( N -V/ 2 det 



a a - bp 



= 1,...,M 



the product term in ea. (|18l) can be represented in the determinant form: 



M 

n 



\\{n a - rip) - i{p a - pp)\ 



M 



a <f3 \K n a +np) -i(p a -Pp)\ 



J\(2Xn a ) det 



Xn a - ip a + Xnp + ipp 



a,P=l,...,M 



(19) 



(20) 



Substituting now the expression for the replica partition function Z(N, A) into the definition of the probability function, 
eg. (ITUt . we can perform summation over N (which would lift the constraint J2 a=1 n a = N) and obtain: 



(-1) 



M 



M=l 



M 

n 



a=l 



dy a dp a 
2vr 



lpAn (s +i) 



n«=l 



det 



-Xn a - ip a + Xnp + ipp. 

(21) 

The above expression in nothing else but the expansion of the Fredholm determinant det(l — K) (see e.g. [24]) with 



the kernel 

K = K[(n,p);(n',p')] 

-oo 



n-lpAn(s+l) 



1 



An — ip + An' + ip' 
Using the exponential representation of this determinant we get 



+00 



dyAi(y+ P ' 2 )(-iy l '^e Xn ' {v+x) 



W(x) = lim cxp|~— — 

M=l 



TrK 



M 



(22) 



(23) 



where 



TrK 



M 



M 

n 



dy a dp 
2tt 



Ai(y a +pl) ^(-lr^e 



n a = l 
1 



(Ani — ipi + An 2 + i^2)(An 2 — ip2 + An 3 + ip%)... (Xum - ipM + Am + ip\) 



Substituting here 



1 



An Q - ip a + An Q+ i + ip a+ i 
one can easily perform the summation over n Q 's. Taking into account that 



du a exp[— (An Q - ip a + Xn a+1 + ip a+1 )uj a ] 



lim V(-l) n - 1 e Anz = lim 



A— J-oo 



i + e 



Az 



and shifting the integration parameters, y a —> y a — x + io a + w ct _ 1 and u> a — > u> a + x/2, we obtain 

M «oo /*+oo j />oo 

lim TrK M = U dy a f2L / du, Q Ai(y Q +j£ + o; a + Wa _ x ) e^""' 



(24) 



(25) 



(26) 



(27) 



where by definition it is assumed that ujq = ujm- Using the Airy function integral representation, and taking into 
account that it satisfies the differential equation, Ai"(i) = t Ai(t), one can easily perform the following integrations: 



oo /"+oo 

dy 



^ Ai(y +p 2 +oj + w'je*^ - "' 5 = 2- 1 / 3 [°° dy Ai(2 1 ^u + y) Aife^V + y) (28) 

27T Jo 

_ Ai(2 1 /3 W ) Ai'(2 1 /3 W ') _ Ai'(2 1 /3 W ) Ai(2 1 /3 W ') 
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Redefining uj a — > u> a 2 1 / 3 we find 



// 



/. 



oo 



lim TrK 



du>icbj 2 ...(kv M K*(loi,u>2)K*(uj2,lo 3 )...K*(lj m ,uji) 



(29) 



-x/2 2 / 3 



where 



Ai(w)AiV) - Ai'(w) Ai(w') 



(30) 



is the so called Airy kernel. This proves that in the thermodynamic limit, L — > 00, the probability function W(x), 
eq.®, is defined by the Fredholm determinant, 



where K* is the integral operator on [—x/2 2 / 3 , 00) with the Airy kernel, cq. (f3T))) . The function ^(s) is the Tracy- 
Widom distribution |9J 



where the function q{t) is the solution of the Panleve II equation, q" = tq + 2q 3 with the boundary condition, 
q(t — > +00) ~ .Az(t). This distribution was originally derived as the probability distribution of the largest eigenvalue 
of a n x n random hermitian matrix in the limit n — >■ 00. At present there are exists an appreciable list of statistical 
systems (which are not always look similar) in which the fluctuations of the quantities which play the role of "energy" 
are described by the same distribution function i^Csh These systems are: the polynuclear growth (PNG) model 
[Tol | . the longest increasing subsequences (LIS) model [ll|, the longest common subsequences (LCS) [13], the oriented 
digital boiling model [l3j . the ballistic decomposition model |l4j . and finally the zero-temperature lattice version 
of the directed polymers with a specific (non-Gaussian) site-disorder distribution Now we can add to this list 
one-dimensional directed polymers with Gaussian <5-correlated random potential. 
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